An existing computer program, used for predicting the natural frequencies and mode shapes of helicopter rotor blades, was refined to improve program accuracy and versatility. The program is based on the Holzer-Myklestad approach adapted for rotating beams. Coupled vertical (out-of-plane), horizontal (in-plane), and torsional mode characteristics were determined for a variety of hub and blade configurations. The resulting program is documented by presenting the recursion equations and techniques for determining natural frequencies and mode shapes, input data requirements, and descriptions of various program outputs. The accuracy of the program is demonstrated by comparing computed results with exact solutions to classical problems and experimental data.